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Abstract. 

■ The Far Ultraviolet Spectroscopic Explorer {FUSE) has surveyed a 

large sample (> 100) of active galactic nuclei in the low redshift universe 
(2: < 1). Its response at short wavelengths makes it possible to mea- 
sure directly the EUV spectral shape of QSOs and Seyfert 1 galaxies at 
z < 0.3. Using archival FUSE spectra, we form a composite extreme ul- 
traviolet (EUV) spectrum of QSOs at 2: < 1 and compare it to UV/optical 
composite spectra of QSOs at higher redshift, particularly the composite 
spectrum from archival Hubble Space Telescope spectra. 



1. Introduction 

The ubiquity with which QSOs display spectral properties such as power-law 
continua and broad emission lines over wide ranges in luminosity and redshift 
has led to the use of composite spectra to study their global properties. Infor- 
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mation about the continuum in the rest-frame ultraviolet is particularly critical 
for understanding the formation of the emission lines, for characterizing the Big 
Blue Bump, and for determining the ionization state of the intergalactic medium 
(IGM). Composite QSO spectra covering the rest-frame ultraviolet have been 
constructed for objects with 0.33 < 2; < 3.6 from HST (Zheng et al. 1997, Telfer 
et al. 2002, T02 hereafter), and at z > 2 from ground-based samples like the 
SDSS (Vandcn Berk et al. 2001), the First Bright Quasar Survey (Brotherton 
et al. 2001) and the Large Bright Quasar Survey (Francis et al. 1991). 

The FUSE bandpass, 905-1187 A, allows us to examine the EUV properties 
of local AGN. We can therefore study the same rest-frame wavelength region 
covered by the HST composite spectra, at rcdshifts less than 0.33. The low 
redshifts of these AGN ensure that, although the FUSE aperture limits it to 
observing relatively bright AGN, our sample contains a large fraction of intrin- 
sically low-luminosity objects. An added advantage to working with low redshift 
spectra is that the determination of the mean EUV spectral index requires a less 
significant correction for IGM absorption than was required for the HST sample. 

2. FUSE Spectra and Composite Construction 

Similarly to T02, we excluded spectra of broad absorption line quasars and 
spectra with S/N ^ 1 over large portions from our FUSE sample. We also 
exclude spectra of objects that show strong narrow emission lines, strong stellar 
features, or strong interstellar molecular hydrogen absorption. A total of 128 
spectra of 90 AGN meet the criteria for inclusion in the sample. We follow the 
same procedure as T02 for the reduction of the sample spectra. To summarize: 
we correct for Galactic extinction using a standard extinction curve, individual 
E{B — V) values for each AGN sightline, and Ry = 3.1; we ignore wavelength 
regions affected by ISM absorption lines; we correct for Lyman limit absorption 
if the S/N below the Lyman break is greater than one; we apply a statistical 
correction for the line of sight absorption due to the Lja forest; we shift the 
AGN spectrum to the rest frame and rcsample to common 1 A bins. 

The lower rcdshifts of the FUSE AGN compared with the HST sample of 
T02 compels us to use different parameters to perform the correction for Lya 
forest absorption mentioned above. Like T02, we describe the distribution of 
absorbers by d'^n/dzdN oc (1 -|- z^N^^. For the column density distribution 
parameter, we use the result found by Dave & Tripp (2001) from echelle spectra 
of two QSOs at z ^ 0.3, (3 = 2.0 for 12.2 < logiV < 14.4. For 14.4 < logA^ < 
16.7, we use /3 = 1.35 from the study by Penton et al. (2000). For the redshift 
distribution parameter, we use 7 = 0.15 (Weymann et al. 1998). We normalize 
the distribution of absorbers by 1.34 x 10~^^ cm^ at logiV = 13 and z = 0.17 
and assume a Doppler parameter of 21 km s~^ (Dave & Tripp 2001). 

We combine the sample spectra using the bootstrap technique described by 
T02. We begin the bootstrap procedure at the central portion of the output 
composite. We then include spectra that fall at longer wavelengths in sorted 
order to longer wavelengths, and spectra at shorter wavelengths in sorted order 
to shorter wavelengths. The overall composite is renormalized at each step. We 
fit a power law of the form F^ oc to the continuum of the composite spectrum; 
and find that the best-fit power law index is a = —0.56. This is equivalent to 
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Figure 1. Top left: Composite AGN spectrum with power law con- 
tinuum fit shown by dashed hue and wavelength regions used in fit 
shown by solid line at bottom; HST composite from T02 shown for 
comparison. Bottom left: Ratio of FUSE to HST composite spectra. 
Right: Composite AGN spectrum with emission line fit. 

the ttEUV fits to wavelengths > 500 A in T02. We show the composite spectrum 
and its ratio to the HST composite in the left panels of Figure 1. On the right, 
we show the emission lines fit to the FUSE composite. 

We find a standard deviation of 0.11 in a from 1000 bootstrap samples of 
the FUSE sample. This gives an estimate of the error arising from the range of 
spectral shapes of the individual AGN that constitute our sample. We explored 
a number of possible systematic errors that could affect the spectral shape of 
the FUSE composite spectrum. The results are sensitive to the extinction cor- 
rection. Changing all individual values of E{B — V) by ztlcr, where we estimate 
la = 0.16E{B - V) (Schlcgcl, Finkbciner, & Davis 1998) changes a by ±0.16. 
Changing Ry in the extinction law to 4.0(2.8) changes a by — 0.19(-|-0.06). The 
composite is also sensitive to the value of the column density distribution pa- 
rameter, p. Reducing it from the chosen fiducial value of 2.0 to 1.7 or 1.5 
increases a by up to 0.3. We estimate the total error from the items above, 
a = — 0.561q28. "^^^ EUV spectral index of the HST composite is significantly 
softer, a = —1.76 ±0.12. In Figure 2, we show the redshift and luminosity 
distributions of the FUSE and HST samples. 

3. Summary 

We summarize our results as follows: 

(1) We construct a composite EUV (630-1155 A) AGN spectrum of objects 
with z < 0.67 from archival FUSE data. (2) We fit a power law continuum and 
Gaussian profiles to the emission lines in the composite spectrum, and wc find 
that O vi/Ly/3 and Ne viii emission are enhanced in the Fi/^^^ composite relative 
to the HST composite due to the Baldwin effect, seen also in the comparison 
between the HST and SDSS composites (T02). (3) We find that the best-fit 
spectral index of the composite is a = — 0.56^o!28- conservative estimate 
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Figure 2. Luminosity versus redshift for FUSE and if5T AGN with 
lines marking redshift and luminosity cuts. 



of the total error in the spectral index includes the standard deviation in a in 
1000 bootstrap samples of the FUSE data set, uncertainties in the extinction 
correction applied to the FUSE spectra, and uncertainties in the column density 
distribution parameter of the intervening Lya forest. (4) The FUSE composite 
is harder than the EUV portion of the HST composite spectrum of T02 who find 
a = -1.76±0.12 for 332 spectra of 184 AGN with z > 0.33. The Baldwin effect is 
generally attributed to the tendency for low-luminosity AGN tend to show harder 
ionizing continua (Zheng &l Malkan 1993; Wang et al. 1998; Dietrich et al. 2002). 
The median luminosity of the AGN in the FUSE sample is logLmedian = 41.2, 
versus log Lmedian = 42.9 for the i7ST sample. One interpretation of these results 
is that both the enhanced high-ionization emission line strengths and the harder 
continuum shape of the FUSE composite spectrum are due to the larger fraction 
of low-luminosity AGN in the FC/S'E' sample. However, we note that splitting the 
FUSE sample itself into high- and low-luminosity/redshift subsamples as shown 
in Figure 2 results in only marginally different values for the EUV spectral index. 
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